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' Starting from the random phase approximation for the weakly coupled multiband 

o 

O ■ tightly-bounded electron systems, we calculate the dielectric matrix in terms of intra- 



band and interband transitions. The advantages of this representation with respect to 
the usual plane- wave decomposition are pointed out. The analysis becomes particularly 
transparent in the long wavelength limit, after performing the multipole expansion of 
bare Coulomb matrix elements. For illustration, the collective modes and the macro- 
scopic dielectric function for a general cubic lattice are derived. It is shown that the 
dielectric instability in conducting narrow band systems proceeds by a common soft- 
ening of one transverse and one longitudinal mode. Furthermore, the self-polarization 
corrections which appear in the macroscopic dielectric function for finite band sys- 



terns, are identified as a combined effect of intra-atomic excliange interactions between 
electrons sitting in different orbitals and a finite inter-atomic tunneling. 
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1 Introduction 

Random phase approximation (RPA) is the well-known textbook method for 
treating the dielectric properties of the weakly interacting electron liquids. Unfortu- 
nately the widely used [2 - 5] extension of this method to the multiband systems leads 
to a cumbersome problem, even when limited only to the calculation of collective modes 
of the tightly bound electrons in the optical regime. In the present text we propose 
therefore a new formulation of RPA which leads to a more convenient description of the 
dielectric response of the tightly bound electrons. In particular it greatly simplifies the 
analysis of the collective modes giving a direct insight into their physical origin. It also 
resolves some long lasting controversies 0, ^ on the local field effects in multiband 
insulators and metals. 

The paper is organized as follows. In Sec. 2 we formulate the problem and derive 
the tight-binding expression for the dielectric matrix. In Sec. 3 the determinant of 
this matrix for the lattice with the cubic symmetry is explicitly calculated in the long 
wavelength limit in which one can make the multipole expansion of the long range 
contributions to the bare Coulomb matrix elements. The collective modes and the 
macroscopic dielectric function which follow from this result are considered in Sees. 4 
and 5 respectively. The latter Section also contains the discussion of the so-called self 
polarization corrections. Concluding remarks are given in Sec. 6. 
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2 Dielectric matrix in the tight binding approach 

In the standard approach p, Q one uses the plane wave representation and starts 
from the infinite dielectric matrix Eq q '(q, uj) with G and G' denoting all the vectors of 
the reciprocal lattice. In order to determine the collective modes and the macroscopic 
dielectric function ej} = (£:^^)oo one has to calculate the determinant of the dielectric 
matrix Eq q >. Unfortunately there was no straightforward way to reduce the dimension 
of this determinant even when only a finite number of presumably relevant tightly 
bound bands was taken into account. The reason is that the plane wave basis used for 
the matrix elements of the Coulomb interaction is not suitable for the description of the 
tightly bound electronic band states. As a consequence, already in the extreme long 
wavelength limit the analysis of the dielectric properties was rather complicated 0, 
while sophisticated and physically nontransparent numerical algorithms were inevitable 
for more complex situations 

One easily appreciates that the tight binding (TB) Wannier-like representation with 
explicit band indices is more appropriate than the plane wave one. Having in mind the 
fact that the basis defined by band indices and the values of wave vector within the first 
Brillouin zone brings in the full crystal symmetry (including the discrete translations) 
like the basis defined in the extended zone scheme, we introduce in the present work the 
dielectric matrix in the former basis, without any reference to the plane wave indices 
G. This approach was already used in studying the dielectric properties of the single 
metallic band 0, but not, to the best of our knowledge, for multiband systems. In 
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fact, the TB basis was frequently used indirectly @, 0, P], i- e. in the calculation 
of the RPA multiband polarizabilities which appear in the matrix elements Bqq >. This 
led to some simplifications in inverting the dielectric matrix Eq q > (q, uj) [|l^ . We point 
out below some crucial advantages of the present direct TB representation. 
We start by calculating the mean values of the components 

pu'{q,t) = 2Y,< a+(k,t)a,,(k+q,t) > (1) 
k 

of the density response to the corresponding components of the external potential 
yiTi^.t) = ^E^'^^/ ^r(r-R)K^^*(r,t)v.r(r-R). (2) 

R, 

Here (/9/(r-R) is the TB orbital for the /-th band at the R-th site, and a/'(k) is the 
creation operator for the corresponding Bloch state. The present discussion is limited 
to the weak coupling in which the bare intraband and interband Coulomb interactions 
are smaller than the bandwidths and the interband energy differences respectively. 
Starting from the equations of motion for the operators of the electron-hole pairs 
a/'(k)aj' (k+q), and performing the standard RPA steps lH, one arrives at the 
system of linear equations 

E(V' - VrMpr = Il,V;^\ (3) 
r 

The indices in the eq.(3) stay for the ordered pairs of band indices, p = {1,1') and 
denote the transitions between the tight-binding bands. Note that the spin indices can 
be omitted in the RPA. The function Ilp{c[,uj) is the RPA polarization diagram 

n,(q,c.) = lv -"'(k+q) (4) 



with the occupation and dispersion of the l-th band given by rii and Ei respectively. 
Evidently, Up is finite only for transitions between (partially) full and (partially) empty 
bands as well as for the transitions within the partially filled bands. 

Vpr{q) in eq.(|^) is the matrix element of the bare Coulomb interaction in which we 
keep only the contributions with two and two TB orbitals centered on same crystal 
sites (two center integrals), 

2 

Vprici) = Y: e^'^^ / d'r [ rfV<(r - R)^^r')^— (r - R)^.. (r'), (5) 
Jrv, 

with p = [li, I'l) and r = {I2, 1'2)- This usual TB approximation which greatly simplifies 
the system of equations (^ and works best for long wavelengths is based on the 
assumption that the overlaps between the orbitals on the neighboring sites are so small 
that the corresponding contributions to the Coulomb matrix elements are negligible 
with respect to the contributions retained in the eq.(^. Basically, this is the only 
restriction on the present RPA approach. In return it leads to simple p, r selection 
rules, as will be discussed below. Further on, in the long wavelength limit it will be 
sufficient to keep only the leading contributions in the multipole expansion of the two- 
site (R 7^ 0) terms in eq.(^. The on-site (R = 0) terms are usually omitted, or at 
least not treated explicitly in the weak coupling limit. However, as it will be argued 
below, some of these terms are relevant for the understanding of the local field effects. 

In the present approach the dielectric matrix is given by the coefficients on the 
left-hand side of eq. (|) , 

^p,r(q, t^) = 5p,r - Vrplip, (6) 

in accordance with the definition which follows from the Dyson equation for the 
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screened Coulomb interaction [jTSl, and with the standard definition for the special 
single band case 

With the wave vector q restricted to the first Brillouin zone and with / covering all 
band indices, the matrix £:p ,.(q, cu) represents, like the matrix Eq q /(q, u;), a complete 
RPA dielectric response, with all local field effects taken into account. Strictly, its 
dimension is infinite, due to the infinite number of possible interband transitions with 
Up 7^ 0. However, an approximation which neglects interband transitions with small 
polarizations Up [e.g. those from the deep (core) orbitals and those to very high 
empty bands], reduces the dimension of ep,r(q, uj) to a finite value, equal to the number 
of retained interband and intraband transitions. This type of simplification with the 
explicit physical justification cannot be achieved in the plane wave representation which 
requires numerical truncations in both band and reciprocal lattice indices ^, 12 . 



3 Long wavelength dielectric response of the cubic lattice 



The further advantage of the dielectric matrix in the TB representation becomes 
apparent in the explicit long wavelength limit. For the sake of definitiveness let us 
consider the cubic lattice and at first take into account only the two-site terms in the 
eq.(^). An analogous procedure can be carried out for any crystal symmetry. We 



distinguish the following three types of dominant lattice sums [|14| in the multipole 
expansion of the various matrix elements Vpr{c[ —>■ 0). 

For both p and r representing intraband transitions [i.e. the transitions within 



partially filled bands, with p = and r = {l2,h)], the dominant contribution is 

the monopole-monopole one, given by |jl5 



VpM) = (7) 

where a is the lattice constant. 

If both indices represent interband transitions [i.e. p = r = (/2) ^2) with 

li 7^ l[ and I2 7^ I2] the lowest possible multipole contribution is the dipole-dipole term 

Here 

t^p = ej dT^^*{r)r^i'{r) (9) 

is the dipole matrix element for the transition from the l-th to the I'-th local function, 
and /ip|| and its respective projections parallelly and perpendicularly to the 

wave vector. 

The third type of matrix elements are those with one interband and one intraband 
transition [i.e. p = {h,l'i) with l[ and r = (hyh) , or vice versa] for which the 
lowest multipole contribution is of the monopole-dipole form, 

_ _«£ ft, ^ 

q 

We emphasize that the long-range part of the Coulomb interaction contributes only 
to the above three types of matrix elements in the long wavelength limit q — > 0. The 
other terms in the multipole expansions are higher orders in components of q, and as 
such are vanishing in this limit. They are therefore irrelevant even if e. g. the dipolar 
factors in the eqs.(§) and ([T0| ) are absent for symmetry reasons. 



Let us now divide all transitions into four sets, taking into account particular prop- 
erties of the cubic symmetry. Each dipolar matrix element (§) is directed along one 
of three fourfold rotation axes (x, y and z) |T^. All interband (dipolar) transitions 
can be therefore divided into three sets, {pj}, where j = 1,2, and 3 stay for x,y and 
z axis respectively. The fourth set {po} contains all (at most three [|T5|) intraband 
(monopolar) transitions. Furthermore, for each dipolar matrix element there are 
two other matrix elements /x^^ and with the same absolute value, which complete 
a three-dimensional irreducible representation of the cubic point group |T^. 

These symmetry properties, together with the factorized form of the matrix elements 
(I) and ([T0|) , reduce the system of eqs.(||) to four linear equations. To this end we 
introduce the following linear combinations of the density components Pp belonging to 
the four sets defined above: 



Po= J2 Pp (11) 

PG{po} 



and 



qn, ^ 

pj = -t — 2^ pppp j = 1,2,3 (12) 

PG{Pj} 



with hj denoting the unit vector along the j—th axis. Summing up separately equations 
(|) with p G {po} and p G {pj},j = 1,2,3, one arrives after few straightforward steps 
to the decoupled systems of two and two equations. The first system represents the 
longitudinal response 



3^2 



[1 + (4vraJ-Vo + Ep, = -^2^'^'^ (13) 



Po + 



2 , . il A . 



| + (47ra,) 1 



EPj=^73-q-Eexi, (14) 
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while two linearly independent equations among the following three equations 



(47ra,)-i - ^ 



quj qrii 

—Pi - — Pj 
qrii qrij 



q X Eext) ■ ilk (15) 



4716 

(with i j ^ k) represent the transverse response. Here q = q/g, y^^'*(q, u;) is 
the external scalar potential and Ee2;t(q, cu) is the external electric field. Eqs.(|T5|) are 
slightly generalized with respect to the original eqs.(3), since they allow for the finite 
transverse electrical field. The intraband and interband polarizabilities are given by 

"c = E n, (16) 



" pe{/} 



respectively. {/} stands for any of three sets {pj},j = 1,2,3. 

The above separation to the longitudinal and transverse response does not depend 
on the direction of q, as it should be for the cubic crystal. Since the original density 
components Pp follow from eqs.(||) once the combinations which figure in eqs.(13 - 15) 
are known, the above procedure gives the explicit solution for the dielectric response 
of cubic crystals. In particular, since the systems of equations (|]) and (13 - 15) are 
connected by linear transformations and thus have a common determinant, it follows 
immediately from the latter equations that 

det[ep^r{cL,uj)] = ei{etf (18) 

with 

ei = {l + Anac) (l + - IGvr^aca/ (19) 
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and 

et = \-^. (20) 

The original determinant of infinite order is reduced so to the algebraic expression in 
which the infinite summation in can be easily truncated due to the explicit physical 
meaning of the terms in eq.(|l^). 

4 Collective modes 



The expressions (18 - 20) are particularly convenient for the discussion of the collec- 
tive modes which appear as isolated poles of det(£) in the cj-plane, and of their Landau 
damping due to the incoherent intraband and interband electron-hole excitations, pro- 
vided the band dispersions and corresponding TB orbitals are specified. The details of 



this calculation for simple band models are discussed elsewhere |^3[. Here we present 
the most important conclusions concerning the collective modes. 

Besides the obvious case of a single band conductor ac 7^ 0, a/ = with intraband 



plasmons as collective excitations, the result (18) covers also a less evident asymptotic 
limit of an atomic (molecular) insulator. Namely, after assuming that all interband 
energies \E\ — Ei\ = Ep are much larger than the corresponding bandwidths, the 
interband polarizability (^7^ reads 

2 UpEplfipl 
pg{/} p 

where Up is the number of electrons per a site that take part in p-th 
transition. Inserting this expression together with = into eqs.( p!9| , pOD , one 
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gets the dipolar collective modes whose spectrum just coincides with that of Frenkel 



excitons [l^] • It should be noted in this respect that the present RPA result is valid in 



m 



the weak coupling regime We remind that the same spectrum is obtained 
the opposite strong coupling regime when the bandwidths are negligible with respect to 
the on-site Coulomb repulsion V^rlR- = 0) between two electrons on different orbitals 
(z. e. the R = contribution in the expression (|^) with li = l[ and I2 = l^)- 

The Frenkel excitonic spectrum thus appears to be a common asymptotic limit of 
the two otherwise incompatible regimes. Indeed, the corrections due to the finite band- 
widths are scaled differently in the two regimes, i.e. by Ep and Vpr{K = 0) respectively. 
In the strong coupling regime the collective modes are coherent superpositions of atomic 
electron-hole excitations, gradually delocalized by finite bandwidths. The appropriate 
representation in the weak coupling regime starts from the electrons in Bloch states. 
Accordingly, in contrast to the former, the latter regime includes also the narrow band 
conductors in which the interband polarizability can be still approximated by eg. (pl|) , 
while the intraband polarizability is finite and given by 

47ra, = -4, (22) 

where Upi = AirUee'^ /m*a^ is the frequency of the intraband plasmon, and m* and rig 
are respectively the effective mass and the number of electrons (or holes) in the metallic 
band(s). 

As can be seen from eq.(|l9D, the finiteness of ujpi causes the renormalization of 



the dipolar longitudinal modes, but does not affect the transverse modes (|20|) . This 



coupling between longitudinal intraband and interband collective modes is to be traced 
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back to the finite monopole-dipole interaction (0). Although the transverse modes (pOD 
are not screened by the intraband charge fluctuations, there is an interesting relation 
between them and the hybridized longitudinal modes of eg. (p!9D . Let the frequency of 



the lowest transverse mode in eg. (pOD , oot^^^, tend to zero. Then the equation ei = 
has the solution, uj, , which also tends to zero. The reverse is also true. The ratio of 

' ''low ' 

two frequencies is 

I ^how \ _ '^pi (23) 



UJt, 2 , 47raj(a;=0) JJ Ep 



where loTj are frequencies of all other (stable) transverse collective modes. Note that 
uji, < UJt, . The dielectric instabilities in multiband conductors are thus characterized 

''low ''low 

by the simultaneous softening of two collective modes. This somewhat surprising, and 
to the best of our knowledge new result is the consequence of the proper treatment of 
the metallic screening of the interband (dipolar) excitations in the present approach. 
We note that the cubic symmetry chosen here is not essential for the validity of this 



result |ll3 |. 



5 Macroscopic dielectric function 



The macroscopic dielectric function em follows immediately from eqs. (p!3|JI^ ). To 
this end it suffices to utilize the averaging scheme by which the probe electrical field 
and the corresponding density response are 'Eiext = iqV"'^'^*(q)/e and p(q) = Z]j=oPi(^) 
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respectively. According to the standard definition Em is then given by 

q, ^ = 1 + 3 ^Zt/A 2 = - = 1 + 4vrae + 24 

This expression reproduces some earher results p, ^ 0, but, unlike them, does not 
contain any type of the so-called self-polarization corrections. In our approach these 
corrections are related to the on-site (intra-atomic) contributions to the bare Coulomb 
matrix elements (|^), not taken into account until now. Obviously, all these contribu- 
tions are independent of q. Furthermore, it is easy to see that due to the symmetry 
reasons only those on-site matrix elements Vpr(R = 0), with p and r representing the 
transitions of the same (monopolar, dipolar, etc) type, are finite. On the other side, in 
the RPA it is sufficient to keep only those on-site terms for which the leading multipole 
contributions in the corresponding two-site sums are constant or vanish. Thus, the 
on-site terms which complete the monopole-monopole two-site sums (|^) in eq(^ [like 
e.g. the on-site "electrostatic" repulsions Vpr(R, = 0) between two electrons sitting 
at the same or different orbitals] are irrelevant in the limit q — > 0. Since the on-site 
contributions which correspond to the monopole-dipole two-site sums ( |T0| ) vanish, the 
most interesting on-site contributions are those going together with the next in order 
( i.e. dipole-dipole) two-site contributions (||). 

The inclusion of those on-site contributions into the expression (H) prevents unfor- 
tunately the reduction of the system (|^) to the explicitly solvable form (^,0). Still, 



some interesting conclusions can be drawn for simple models [0, like that with only 
two bands (/ = hjh) connected by a finite dipole matrix element (j^) fi.p,p = {h,h)- 
Then the dipole-dipole matrix element Vp^p of eq(||) has to be completed by the on-site 
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term 



V,,piR = Q) = jdhj rf3/y,;^(r)y,^r')^:^¥.Jr)y,,,(r') (25) 

which has the meaning of the exchange between two atomic orbitals. Furthermore, the 
interband polarizabihty ( p!7D reduces to 

aI = -\\^lWIlp + Ilp). (26) 

The macroscopic dielectric function Em can be still written in the form (0), but with 
the interband polarizabihty a/ replaced by 

I Oil 

= i-i;,^(i? = o)(n, + n^) • ^^^^ 

Note that there is no effect of this type on the intraband (metallic) polarizabihty etc, in 
accordance with the already mentioned fact that for small q the long-range monopole- 
monopole Coulomb interaction dominates over all on-site "electrostatic" repulsions. 

The denominator in the eq.(^) can be simply interpreted as the atomic screening 
of the interband polarizabihty a/. Since it has itself the RPA form, the formulation 
(i) apparently treats both, the band and the local (intra-atomic) dielectric responses 
selfconsistently at the same (RPA) level of approximation. 

When the bandwidths are finite, the intra-atomic and the inter-band screenings in 
eg . (|27D cannot be decoupled. As a consequence the result for Em deviates from the 
Lorentz-Lorenz (LL) form in the sense that the effective interband polarizabihty olj 
which enters into eq (^4]) cannot be reduced to the form (p!7|), i. e. to the sum of 
interband polarization diagrams. This deviation is usually named a self-polarization 
correction, originally derived by Adler within the plane-wave representation. The 
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result (^) brings us to the clear microscopic interpretation of this correction, recog- 
nized as a combined effect of the intra-atomic screening due to the on-site exchange 
Coulomb interactions and of the interatomic tunneling {i. e. of the finite bandwidths). 
Note furthermore that some later results for the macroscopic dielectric function P, 
cannot be represented in the form (pTI). The reason might be in the additional approx- 
imations of these works which, in contrast to the present approach, do not include fully 
the local RPA on-site screening into the dielectric matrix. 

In the asymptotic limit of zero bandwidths the expression for sm reduces again 
to the LL form, but with the initial (e. g. Hartree-Fock) values of the dipole matrix 
element /i and the energy difference Ep replaced by |/ie//P = |/^P-E'p/-E'p,e// and E^ eff = 
Ep + 2neEpVp^p{R = 0) respectively. In the present approach this renormalization is 
performed within the RPA scheme, i. e. the result (p7[) includes both the crystal and 
the on-site dielectric screening on the same RPA level. It is appropriate to remind 
here that with vanishing bandwidths such approach is valid only if the interactions 
are small with respect to Ep. It can be also noted that Adler uses the term "self- 
polarization correction" for the on-site contributions to the dielectric screening even in 
the zero-bandwidth limit. 

The inclusion of the intra-atomic Coulomb interactions does not alter the previous 
result (|23| ) concerning the common instability of the transverse and longitudinal modes 
in narrow band conductors. E. g, in the asymptotic limit of zero bandwidths we get, 
after the inclusion of the on-site "exchange" Coulomb matrix elements into the RPA 
calculation, both the collective excitations and the macroscopic dielectric function of 
an atomic insulator with the RPA values of the intra-atomic parameters. If the latter 
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represent an appropriate description of the on-site correlations, our method is aheady 
fully self-consistent. In the case of the strong on-site correlations, the short-range 
contributions to the dielectric response have to be treated beyond RPA. 

Finally, in the limit of large frequencies (a; — oo) the macroscopic dielectric function 
(p^) reduces to the simple sum of intra-band and inter-band polarizabilities, with the 
standard free electron mass expression for the plasma edge 

6 Conclusion 

In conclusion, we point out the main advantages of the present TB approach to 
the dielectric response. Unlike the representation via the reciprocal lattice indices, 
it leads to a simple and physically transparent expression for the dielectric matrix, 
which is explicitly resolved in the particular example of the cubic symmetry. It also 
clearly distinguishes between the long-range and the local (on-site) contributions to the 
dielectric screening, and gives so the direct insight into the origin of the self-polarization 
effects. Furthermore, the present method not only interpolates between (multiband) 
conductors and insulators, but also includes systems with neghgible bandwidths as a 
well-defined asymptotic limit. Just in this limit we discover an important property 
of the dielectric instability in multiband conductors, namely that it proceeds by a 
simultaneous softening of one longitudinal and one transverse collective mode. 

In summary, the approach proposed here facilitates the detailed analysis of the 
dielectric response in crystals. It is particularly efficient in the long wavelength limit, 
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even if a large number of bands has to be retained in the calculations. In this respect 
the decisive criterion for the reducibility of the problem to a finite number of linear 
equations is the smallness of the on-site exchange terms (^) in comparison with the 
corresponding long range contributions. Other on-site terms are not relevant (but they 
have to be weak enough in order to justify the starting RPA scheme). On the other 
hand, it is clear that our approach has a wide range of applications to real and model 
systems in which it is appropriate to take into account a small number of bands. 
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